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Abstract 

We present a simple time dependent model for the excitation of 
a nucleon from a bound state to a continuum resonant state in a 
neutron-core complex potential which acts as a final state interaction. 
The final state is described by an optical model S-matrix so that both 
resonant and non resonant states of any continuum energy can be 
studied as well as deeply bound initial states. It is shown that, due to 
the coupling between the initial and final states, the neutron-core free 
particle phase shifts are modified, in the exit channel, by an additional 
phase. The effect of the additional phase on the breakup spectra is 
clarified. As an example the population of the low energy resonances 
of ^^Be and of the unbound ^^Be is discussed. Finally, we suggest that 
the excitation energy spectra of an unbound nucleus might refiect the 
structure of the parent nucleus from whose fragmentation they are 
obtained. 
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1 Introduction 



In this paper we will call projectile fragmentation the well known elastic 
breakup (diffraction reaction) of neutron halo nuclei, when the observable 
studied is the neutron-core relative energy spectrum. This kind of observable 
has been widely measured in relation to the Coulomb breakup on heavy 
target. Recently results on light targets have also been presented [T]. These 
data enlighten the effect of the neutron final state interaction with the core 
of origin, while observables like the core energy or momentum distributions 
enlighten the effect of the neutron final state interaction with the target. 

Projectile fragmentation has been used experimentally also with two neu- 
tron halo projectiles. In this case it has been suggested that the reaction 
might proceed in one step (simultaneous emission of the two neutrons) or 
two steps (successive emissions) depending on whether the target is heavy 
and therefore Coulomb breakup (core recoil) is the dominant mechanism or 
the target is light and then nuclear breakup is the dominant mechanism [2]. 
The successive emission can be due to different mechanisms. One possibility 
is that one neutron is ejected because of the interaction with the target, as 
in the one-neutron fragmentation case, while the other is left behind, for 
example in a resonance state, which then decays. This second step has been 
described by the sudden approximation in Ref. j3| under the hypothesis that 
the first neutron is stripped and that the transparent limit for the second 
neutron applies. It corresponds to consider the second neutron emitted at 
large impact parameters such that the neutron-target interaction can be ne- 
glected. The two-step mechanism implies that the two neutrons are not 
strongly correlated such that the emission can be considered sequential. 

However the neutron-target interaction gives rise not only to stripping 
but also to elastic breakup and in both cases to first order in the interaction 
the neutron ends- up in a plane wave final state It can then re-interact 
with the core which, for example, is going to be ^°Be in the case of the one- 
neutron halo projectile ^^Be, while it will be ^^Be in the case of the projectile 
fragmentation of ^"^Be, since ^^Be is not bound. While in the case of ^^Be the 
structure of both its bound and continuum states is well known from other 
kinds of experiments and therefore projectile fragmentation experiments are 
useful to enlighten the reaction mechanism and its possible description, in 
the case of ^^Be or of other unbound nuclei the interplay between structure 
and reaction aspects is still to be clarified. 

Experiments with a ^^B projectile [3 E] have also been performed, in 
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which the n-^^Be relative energy spectra have been reconstructed by coinci- 
dence measurements. In such a nucleus the valence neutron is weakly bound, 
with separation energy S„=0.969 MeV, while the valence proton is strongly 
bound with separation energy 8^=17.3 MeV. Thus the neutron will proba- 
bly be emitted in the first step and then re-scattered by the core minus one 
proton nucleus. The projectile-target distances at which this kind of mecha- 
nism would be relevant are probably not so large to neglect the effect of the 
neutron-target interaction. One might wonder therefore on how to describe 
a neutron which breaks up because of the interaction with the target, is left 
in a plane wave moving with the same velocity of its original core and re- 
interacts with it in the final state. This mechanism could be at the origin 
of the coincidence measurements for a one-neutron halo system like ^^Be or 
for a projectile like ^^B. It could be also one of the mechanisms giving rise 
to ^^Be in fragmentation measurements of ^^Be. Supposing the two neutrons 
strongly correlated and being emitted simultaneously due to the interaction 
with the target, the coincidence measurement of one neutron with the core 
would evidence the neutron-core final state interaction. 

Light unbound nuclei have attracted much attention [T]-[2B] in connection 
with exotic halo nuclei. Besides, a precise understanding of unbound nuclei 
is essential to determine the position of the driplines in the nuclear mass 
chart. In two-neutron halo nuclei such as ^He, ^^Li, ^'^Be, the two neutron 
pair is bound, although weakly, due to the neutron-neutron pairing force, 
while each single extra neutron is unbound in the field of the core. In a 
three-body model these nuclei are described as a core plus two neutrons. The 
properties of core plus one neutron system are essential and structure models 
rely on the knowledge of angular momentum and parity as well as energies 
and corresponding neutron-core effective potential, therefore spectroscopic 
strength for neutron resonances in the field of the core. Ideally one would 
like to study the neutron elastic scattering at very low energies on the " core" 
nuclei. This is however not feasible at the moment as many such cores, like 
^Li, ^^Be or ^^B are themselves unstable and therefore they cannot be used 
as targets. Other indirect methods instead have been used so far, mainly 
aiming at the determination of the energy and angular momentum of the 
continuum states. 

Unbound nuclei have been created in several different ways besides the 
projectile fragmentation [2l|5l[6l|8]-[T7j mentioned above: multiparticle trans- 
fer reactions [18]-[23] or just one proton [25] stripping. In a few other 
cases the neutron transfer from a deuteron [2S]-[2H] has been induced and 
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the neutron has undergone a final state interaction with the projectile of, 
for example ^^Be. In this way the ^^Be resonances have been populated in 
what can be defined a "transfer to the continuum reaction" [29]-|33]. Thus 
the neutron-core interaction could be determined in a way which is somehow 
close in spirit to the determination of the optical potential from the elas- 
tic scattering on normal nuclei. In both the projectile fragmentation or the 
transfer method the neutron-core interaction that one is trying to determine 
appears in the reaction as a "final state" interaction and therefore reliable 
information on its form and on the values of its parameters can be extracted 
only if the primary reaction is well under control from the point of view of 
the reaction theory. 

In a recent paper [7] we showed that among the methods discussed above 
to perform spectroscopy in the continuum, the neutron transfer method looks 
very promising since the reaction theory exists and it has been already tested 
in many cases [7], [29]- [33]. It is important to remember that the final state 
interaction of the neutron with the target (or with the projectile, in the case 
of inverse kinematics reactions) is contained in the transfer to the continuum 
method developed in Refs.[25]-[55]. 

In this paper and in particular in Sec. 2 the basic formalism to describe 
projectile fragmentation, an inelastic-like excitation to the neutron-core con- 
tinuum [311 ES], is presented and the effect of final state interaction of the 
neutron with the projectile core is studied. The model is a theory which 
would then be relevant to the interpretation of neutron-core coincidence mea- 
surements in nuclear elastic breakup reactions. In the present work we apply 
it to the breakup of the halo nuclei ^^Be, ^^B but also ^^Be. In the case of 
two nucleon breakup we try to describe here only the step in which a neu- 
tron is knocked out from the projectile by the neutron-target interaction to 
first order and then re-interacts in the final state with the core. The case 
in which a resonance is populated by a sudden process while the other neu- 
tron is stripped has been already discussed in Ref.[3| and we will show that 
there is a simple link with the model presented here. The present model is 
therefore partially related to Ref.[3]. We assume that the neutron which is 
not detected has been stripped while the other suffers an elastic scattering 
on the target. But while in Ref. [3] the so-called transparent limit was used 
for the second neutron, corresponding to no interaction at all between the 
neutron and the target, we will consider here explicitly the effect of such 
an interaction on the n-core relative energy spectrum. This will result into 
a core-target impact parameter dependence for the fragmentation form fac- 
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tor. However, in most of our calculations, we shall also use the no-recoil 
approximation for the core (cf. Fig.l). On the other hand the influence of 
a possible second nucleon, when appropriate, is taken into account only by 
a modification of the neutron-core interaction in the final state. A simple 
idea for relating the present work to its future development into a two nu- 
cleon breakup model is presented in Sec. 3. Section 4 contains the results 
of our numerical calculations for ^^Be which, being already well understood, 
has been used here as a test case. It also summarizes experimental results 
and the present theoretical understanding of ^^Be. Furthermore details on 
our assumptions for the potentials needed in the calculations are presented. 
Numerical results for ^^Be are contained in Sec. 5. Finally our conclusions 
are contained in Sec. 6. 

2 Inelastic excitation to the continuum. 




Figure 1: Coordinate system used in the calculations 

To first order the inelastic-like excitations can be described by the time 
dependent perturbation amplitude [HEHES]: 

1 c°° 

Af. = - dt{^Pf{r,t)\V,{r-K{t))Mr,t)), (1) 
in J —oo 

for a transition from a nucleon bound state ipi to a final state ipf which 
can be a bound state or a continuum state. In this paper we shall treat 
only continuum final states. V2 is the interaction responsible for the neutron 
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transition (cf. Eq. (2.15) of [3]). The potential V2(r — R.(t)) moves past on 
a constant velocity path with velocity v in the z-direction with an impact 
parameter he in the x-direction in the plane y = 0. These assumptions 
and the other discussed in Sec. 2.4 make our semiclassical model valid at 
beam energies well above the Coulomb barrier. This is in fact the regime in 
which projectile fragmentation experiments are usually performed (cf. Sec. 
4). The coordinate system used in the calculations is shown in FigJT] and 
it corresponds to the no-recoil approximation for the core. In the case of 
the very weakly bound ^^Be we will drop this approximation and explicitly 
take into account core recoil by defining R(t) as the projectile-target relative 
motion coordinate. 

Let ipiir^ t) = 0j(r)e~fi'^'* be the single particle initial state wave function. 
Its radial part 0i(r) is calculated in a potential Vws{^) (cf- Sec. 4) which is 
fixed in space. In the special case of exotic nuclei the traditional approach 
to inelastic excitations needs to be modified. For example the final state 
can be eigenstate of a potential Vi modified with respect to Vws because 
some other particle is emitted during the reaction process as discussed in the 
introduction. The final state interaction might also have an imaginary part 
which would take into account the coupling between a continuum state and 
an excited core. For these reasons our treatment is closer to the formalism 
used in Ref. |35] to treat rearrangement collisions between heavy ions than 
to the inelastic excitation formalism of [51]. For initial and final states of 
different angular momentum our wave functions are trivially orthogonal due 
to the orthonormality of their angular parts. For transitions conserving the 
angular momentum of the single particle states we can orthogonalise the 
initial and final states as suggested by Eq. (8), pag. 303, Ch.V.3, of Ref. 
[Sg. Then 

and the first order time dependent perturbation amplitude reads 

Af,, = - / dtdr 0}(r)0,(r)e*-V2(r - R(t)) - AAj,, (3) 

III J— oo 

where hut is the energy difference between the initial and final states and 

AAfi = -(0/10.) / dtdr 0*(r)0,(r)e*-*r2(r - R(t)). (4) 

til J —oo 

Because of the displacement of V2 with respect to 0j and in particular for 
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the choice of a 5-potential discussed below, it is expected that the correction 
due to AAfi will be small. Therefore we shall neglect it in the following. 

Now change variables and put z — vt = z'oYt = {z — z')/v. The excitation 
amplitude becomes 



^/i = dxdydzdz' (f)*f{x, y, z)(f)i{x, y, z)e"^^^ "-'^V^ix - y, z'), (5) 



where 



Then 



1 

Afi = — dxdydz (t)){x, y, z)(j)i{x, y, z)e"^''V-2{x - be, y, q), (7) 

liXV J — oo 

where 

/oo 
dzV2{x-be,y,z)e'i^ (8) 
-oo 

In our approach the presence of the target represented by this interaction has 
the effect of perturbing the initial bound state wave function and allow the 
transition to the continuum by transferring some momentum to the neutron. 
For this purpose, although the potential V2{r) has a radius of the order 
of the potential of the target, it is enough to choose a simplified form of 
the interaction. Therefore we choose V2(r) to be a delta-function potential 
V2(r) = V2S{x)5{y)6{z), with V2 = [MeV fm^]. Then the integrals over x and 
y can be calculated giving 

^/^ r '^K^c, 0, z)Mbc, 0, z)e"i\ (9) 

tllV J —oo 

The value of the strength f 2 used in the calculation is discussed in Sec. 5 and 
in Appendix A. From the above equation it is clear what the effect of the 
n-target ^-interaction is in a time dependent approach: while in the sudden 
approach the initial and final state overlap is taken in the whole coordinate 
space, irrespective of the target and of the beam velocity, here the overlap 
of the initial and final wave functions is taken at the core-target impact 
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parameter distance on the x-direction which is along the distance of closest 
approach. The y component is zero (neutron emitted on the reaction plane 
preferentially) and the z-component, being along the relative velocity axis is 
boosted by a momentum q. 

The delta-function potential should be a good approximation if 7RT << 1 
and kRx « 1 where Rt is the radius of the target, 7 = -^—2mei/h is the 
decay length of the initial single particle wave function corresponding to the 
nucleon binding energy Si and k = yJ'lmSf/Ti is the nucleon final momentum 
in the continuum state. In this case the initial and final wave functions are 
rather constant over the volume of the target and can be replaced by their 
values at the center. The second condition is related to the first because the 
cross section becomes small if k is large compared to 7 (see Sec. 2.3). We also 
require that the reaction should be peripheral in the sense that Rc + Rt < be 
where Rc is the projectile core radius. The first condition 7Rr << 1 means 
that the projectile should be a good halo nucleus. Another situation where 
the delta-function potential is a good approximation is when 7RT >> 1. 
That is the initial state should be strongly bound and the initial state wave 
function decays rapidly inside the target. Then the (5-potential should be 
located at the surface of the target. If, for example, V2(r) is a square well 
potential with radius Ry and if 7Rt' >> 1, then Eq.([9]) can still be used to 
estimate the breakup with the following changes: 

• be is replaced by be — Rt i-e. the interaction is located at the surface 
of the target. 

• V2 is replaced by V2 where 



Thus the strength is reduced and there is an extra phase. The derivation of 
Eq. ( ITOi) is presented in Appendix A. 

First we study the simple case where the initial bound state and the final 
continuum state have k = If = 0, then 



3 

V2 = - 



2(7-ifc)2i?2 • 



(10) 



Ci e-T"" 



(11) 




r 




Cf k 



:i-{h}^'\kr) - Shl)^>{kr)). 



(12) 
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These expressions are the asymptotic forms of the initial and final state wave 
functions. Their use can be justified when the impact parameter is sufficiently 



large [3j and r = yb^ + z"^. 7 and k are the neutron momenta in the initial 
and final states already defined. Cj is the asymptotic normalization constant 



of the initial state wave function while Cf = \j'^/L is the normalization 
constant for the final state. L is a large box radius used to normalize the 
continuum wave function (cf. Eq. (2.5) of Ref. [29]). The quantity S is 
the S-matrix representing the final state interaction of the neutron with the 
projectile core. Then 

Afi = -- ^ — / dz cosg^;. (13) 



Let us define 



and 



such that: 



while 



then 



and 



Where now C 



roD g-(7-ifc)r 

lR = Re / dz cosqz, (14) 



roo Q-{y-ik)r 

Ij = Im / dz cosg^;, (15) 

J-00 



I{k,q) = lR + tIj = \I\e"' (16) 



S = Se^'" = e^'^^^"^ (17) 



Afi = C{I - S*r) (18) 



Af,\^ = C^\I\^\1- S\^. (19) 



hv 8lT 
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2.1 Wave functions for general / 



In the general case of a Zj > initial state the amplitude Eq.([T]) has to be 
calculated numerically using, for example, the following forms for the wave 
functions. For the initial bound state 

0,(r) = -azS/^l'^(^7r)l/„„.,(^^,0). (20) 

Due of the strong core absorption discussed in Sec. 2.4 and to get a simple 
insight at the physics of unbound nuclei, we use in this paper the asymptotic 
form of the initial state wave function, however the exact wave function, nu- 
merical solution of the bound state Schrodinger equation can be used without 
introducing further complexity in the calculations. 
For the final continuum state 

<Pf{r) = Cfk'-{h\-\kr) - Si,.h[;\kr))Yi^.,^^.{e,</^). (21) 

As it was shown in Ref.[29], in the case of narrow isolated resonances the 
treatment of the continuum states via the S-matrix is equivalent to the R- 
matrix formalism. 



2.2 Probability spectrum 

The probability to excite a final continuum state of energy is an average 
over the initial state 

^in = 21+1 ^™''™/ 

and a sum over the final states. Introducing the quantization condition 

kL = nn, (23) 
and the density of final states, according to Ref. [29] 

pief)def = ^^def, (24) 

TT h k 

the probability spectrum reads 

dPjn _ 2 f I 2 ^ 1 V h _ Q l^lr |2 

dSf 7rZlV^/l2fc2/, + l "-"^^l ^m,,m,\ \lm,,mf\ , [^^ ) 
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where now 



\-^mi,mf I 



■ (26) 



For simplicity the equations in this section are obtained without spin 
variables in the initial and final states. The generalization including spin is 
given in Appendix B. 



2.3 Approximate evaluation of the integral I{k, q) 

In order to study the qualitative effects of the final state interaction we 
proceed now to an approximate evaluation of the integral I{k, q) for /j = 
If = 0. However the calculations presented in Sec. 5 use the exact integrals. 
For large impact parameters be, write r = ^b'^ + z"^ ^ 6^, + /2bc- Then 

0^ J — oo 

Hence the phase u will be given by 



1 f -u^^uh ^^^^^ 
u = — CLrgy") — tk) + kbc — 



2 ' 2(72 + F) 

= --arg{^-tk) + kb,^l-l^j (28) 

where k = mv/h and we have used Eq.([n]) to obtain q = (7^ + k'^)/2k. The 
estimated value of |/P is 



1 27r6c _27fe, 

=e ' " exp 




' ' &^ yJi' + k'^ 

The above analytical expressions are accurate to within 10% for impact pa- 
rameters around the strong absorption radius and for neutron-core energies 
less that l.SMeV. The agreement improves for larger impact parameters. 
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The approximate formulae give rise to simple physical interpretations. 
The first is that we have an explicit expression for the dependence of v on 
the neutron wave number fc, the core-neutron impact parameter he and the 
adiabaticity parameter qjk. We will discuss some of the effects of the new 
phase V in Sec. 5. |Jp ~ can be interpreted as an inelastic-like form 

factor and it is interesting to compare it to the transfer to the continuum 
form factor given in Ref.[29]. The inelastic form factor decreases with 

the impact parameter much faster than the transfer form factor. This is a 
well known characteristic for final bound states [35] and it is interesting to see 
that it persists for final continuum states. Furthermore the slope parameters 
are in both cases given in lowest order by the initial state decay length 7. 

Finally we make connection with the sudden approximation formula Eq.(20) 
of Ref . [3] which describes the second step of a two neutron breakup reaction 
as a resonance decay, when the first neutron has been stripped. In our nota- 
tion it reads 



do \ (k cos 5 + 7 sin 5 \ ^ 

&7 ~ A;(72 + F) 1^ + ) 

1 |sin(5 + /?)|2 

k {'j^ + k'^) 



(30) 



where f5 = arctan{k/'~f). This formalism also predicts the presence of an 
extra phase shift (3 with respect to the free particle scattering determined by 
S since | sin(5 + = ||1 — S\'^ and S = e^^^^^^\ Similarly to our case the 
effect of P would be to modify the resonance-like structures. In both cases 
then S and S could be interpreted as off-the-energy-shell S-matrices. On the 
other hand our additional phase u contains an explicit dependence on the 
impact parameter and we calculate the potential phase shift and S-matrix 
by an optical model code. Our S-matrix can in principle be complex to allow 
for core excitation effects. Also it can consistently, and in the same formal- 
ism, describe resonant and non-resonant final continuum states of angular 
momentum If = but also If > 0. In the latter case, the Breit-Wigner 
assumption for the line shape of the resonances used in other approaches, 
is naturally given by the optical model calculation of the factor |1 — Sl"^ in 
Eq.(l25D. 
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2.4 Cross section 




In [35] it was shown that the semiclassical treatment of peripheral quasi- 

elastic reactions is valid for transfer reactions as well as for (inelastic) pro- 
jectile excitation and therefore we will apply it in the following by simply 
substituting Pin{bc) to Pt{bc) in the well known formula which gives the cross 
section in terms of the neutron excitation probability and the core elastic 
scattering probability. A full description of the treatment of the scattering 
equation for a nucleus which decays by single neutron breakup following its 
interaction with another nucleus, can also be found in Refs. [U [291 ESI [38] 
where the cross section differential in Sf, the final, continuum, neutron energy 
is given as 



(see Eq. (2.3) of [33]) and C^S is the spectroscopic factor for the initial state. 

The core survival probability Pctipc) = \Sct\'^ [33] in Eq. f[3T|) takes into 
account the peripheral nature of the reaction and naturally excludes the 
possibility of large overlaps between projectile and target. Pd is defined in 
terms of a S-matrix function of the core-target distance of closest approach 
be- A simple parameterisation is Pct{bc) = ^^'^^^^^^^"'^"''^"'^^ [33J, where the 
strong absorption radius ~ lA{Ap^^ + A]^^) fm is defined as the distance 
of closest approach for a trajectory that is 50% absorbed from the elastic 
channel and a=0.6 fm is a diffuseness parameter. The values of R^ thus 
obtained agree within a few percent with those of the Kox parameterization 



Because Pinipc) depends on b^ trough the form factor Imi,mf, the final 
cross sections will get the main contribution from a limited range of impact 
parameters around the strong absorption radius. As we shall see in the 
following (cf. figs. 7 and 8) this localization makes the tails of the spectra 
from Eqs. fl25l) and (I3T1) to decay faster than Eq. (!30l) towards neutron-core 
high energies. 

3 Two-nucleon breakup 

This section introduces the basis for a generalization of Eq.f l3ip which will be 
further expanded elsewhere. The goal is to take into account explicitly the 




(31) 



m 
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presence of a second, stripped, nucleon as in the cases of projectiles like ^ Be 
which have a two-neutron halo or ^^B in which a proton can also emitted. 

We start from hypothesis similar to those leading to Eq.(7) of Ref. |3] 
which gives, in the case of two-neutron breakup, the one-neutron-core rela- 
tive momentum distribution when the other neutron is stripped. Let us call 
(1) the stripped neutron and (2) the neutron detected in coincidence with the 
core. Following the approximations proposed by Ref.[3j for the coordinate 
variables, shown in Fig|2l of the core and neutrons with respect to the target 
one gets: Ri±=ri_L+bc and R2±=r2±+bc, where the heavy-core approxima- 
tion has been used, ri and 12 are the coordinates of neutron (1) and (2) with 
respect to the core, while Ri and R2 are with respect to the target. 







































> 






X 





Figure 2: Coordinate system used in the calculations of two-nucleon breakup. 

Suppose the initial two-neutron wave function to be given by a prod- 
uct of single particle wave-functions, as in the shell model: \E'(ri,r2) = 
ai[0i(ri)02(r2)] with spectroscopic factor of = C'^S. For simphcity we con- 
sider here only a = = transition and thus we do not include spin wave 
functions. Then the neutron-core cross section differential in the relative 
energy is: 

^ = 2C'S I d%^\S^,{W)\'^^^ I dh,^{l-\Sr{K,^r) I dz\cl^,{v,^,z)\\ 
CLEf J asf J J 

(32) 

Where dPin{bc)/dsf is given by Eq. (!25|) and, similarly to Eq. (l3T]) . we are 
treating the core-target interaction in the eikonal approximation. The d'^ri± 
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integral above gives the neutron (1) stripping probability. For the Si-matrix 
we consider a sharp cut-off approximation such that S'i(Ri^) = if < 
Ri_L < Rt, while 5'i(Ri_l) = 1 if Rij_ > R-r and Rt is the target radius. 

Thus we obtain a simple expression for the two nucleon breakup cross 
section, in which one is stripped by the target while the other is elastically 
scattered and interacts with the core in the final state 



^ = 2C'S I d%.\S^M\'^^^ I dh.^l dz\Mri±,z)\^ (33) 

as a product of the neutron (2)-core relative energy distribution and a factor 
depending on the stripped neutron (1) wave function. For each core-target 
impact parameter be the limits of the integral on ri_L are: be — Rt < fi± < be. 



4 Applications 

4.1 The reaction ^^Be n+i^Be 

As a test of our model we calculate the relative energy spectrum n+^°Be 
obtained by the authors of Ref. [1] in the breakup reaction of ^^Be on ^^C at 70 
A.MeV. The structure of ^^Be is well known: the valence neutron is bound by 
0.503 MeV; the wave function is mainly a 2s state with a spectroscopic factor 
around 0.8 and there is also a small 6.5/2 component. The main 6.^/2 strength 
is in the continuum centered around 1.25 MeV [37]. We have calculated the 
initial wave function for the s-state in a simple Woods-Saxon potential with 
strength fitted to the experimental separation energy and whose parameters 
are: ro=1.25 fm, a=0.8 fm. As possible final states we have considered 
only the s, p and d partial waves calculated in the /-dependent potentials of 
Table [H The delta-function potential strength has been chosen as -4057.59 
MeV fm^. The authors of Ref.[l] have shown that the effect of Coulomb 
breakup is noticeable in their n+^°Be spectrum. We have also included this 
contribution, calculating it with the method of Margueron et al. [3H] which 
explicitly takes into account the recoil of the core. 

The spectrum of Figl3]is very similar to the spectrum obtained in Ref. [39] 
by solving the time-dependent Schrodinger equation numerically, expanding 
the projectile wave function upon a three-dimensional spherical mesh. Sim- 
ilarly to the present model, a classical, straight line trajectory for the core- 
target scattering was used in Ref. |39j . Also our n-core potentials are very 
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Figure 3: n-^°Be relative energy spectrum, including Coulomb and nuclear breakup 
for the reaction ^^Be+^^C n+^"Be+X at 69 A.MeV. Only the contributions from an 
s initial state with spectroscopic factor C^S= 0.84 are calculated. The triangles are the 
total calculated result after convolution with the experimental resolution function. The 
dots are the experimental points from Ij. 



close to those of Ref.[39] and our ^-interaction strength is consistent with the 
volume integral of their neutron-target interaction. 

We have then folded the calculated spectrum througth the experimental 
resolution function of Fukuda et al. [1], as given in Ref.[39] 



de J dSf 



g{ef-e)def 



with Ef the energy in the theoretical calculation. The result is shown in FigO] 
by the triangles. The full curve is the total spectrum, sum of Coulomb and 
nuclear breakup. Each individual transition, due to the nuclear interaction 
only, is also shown. The dots are the experimental points from [Ij. Our 
calculations include an initial state spectroscopic factor C^S=0.84. The kind 
of discrepancy between our calculation and the data in the range 1-2 MeV is 
very similar to that of the calculations in Ref . |39j . 
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Table 1: Woods-Saxon potential parameters for the s, p and d states in the ^^Be contin- 
uum. 





/=0 1=1 1=2 


V,o(MeV) 


a (fm) 


R (fm) 


Vo(MeV) 


-62.52 -39.74 -63.57 


5.25 


0.6 


2.585 



Encouraged by the good agreement of our results with those of Refs.[T] 
and [39J, we conclude that our model is quite reliable to calculate projectile 
fragmentation including final state interaction with the core and we turn now 
to the study of more challenging reactions. 

4.2 The reactions ^^Be^^^Be +n and ^"^B^^^Be +p 

An interesting unbound nucleus which has recently attracted much attention 
is ^^Be. It has been obtained in several different type of experiments with 
normal and exotic beams but its structure is not clear yet. 

One of the aims of this paper is to see whether the neutron-^^Be relative 
energy spectra obtained from fragmentation of ^"^Be or ^^B would show dif- 
ferences predictable in a theoretical model. Therefore we start by describing 
briefly the present knowledge of ^^Be. 

4.3 Structure of ^^Be 

The first experimental evidence of ^^Be was recorded in Ref.|18]. It was 
the unobserved particle in the two-body reaction ^^C(^Li,'^B) at 82 MeV. 
A narrow resonance at 2 MeV above the neutron emission threshold was 
observed in the spectrum of the measured '^B ions and it was interpreted as 
being due to the ground state of ^^Be. Later on it was identified as a 6.^/2 state 
in the double-charge exchange experiment ^^C(^'^C/^0) at Ej„c=337 MeV 
|22j and in an inverse kinematics (d,p) transfer reaction at 55 A. MeV [26] : 
both the missing mass method, from the detected proton spectrum, and the 
invariant mass spectroscopy, from the measurement of all the decay products 
from the unbound state were used. The 6.^/2 resonance was considered as the 
ground state of ^^Be until the experiment of Ref. [23] • This experiment used a 
stable beam multi-nucleon transfer process ^^C(^^B,^^N)^^Be in which a lower 
state, unbound by 800 keV was observed. A spin J=l/2 was suggested but 
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without parity assignment. More recently a broad peak has been obtained 
in several projectile fragmentation experiments [51 EllISl US] and tentatively 
identified as a 1/2+ state. The experiment of Ref. p0] used fragmentation of 
^^O at 80 A.MeV. Neutrons in coincidence with ^^Be were detected and the 
broad peak was observed in their relative velocity spectrum. The spectrum 
was fitted with a virtual s-state of energy 60 keV while in Ref.[S], which used 
fragmentation of ^^B a virtual s-state could not fit the experimental neutron 
spectrum. The assumption of an s-resonance of energy 700keV and width 
1.3 MeV leads to a good agreement. In Refs.[5l El [151 [IS] these unbound 
states of ^^Be have been populated by breakup of ^^Be or ^'^B. Both types of 
experiments show a low energy peak but only in the breakup of ^^B [51 [S] the 
other peak at 2 MeV corresponding to the 6.^/2 state, is seen clearly. Finally 
Simon et al. [I5| have fitted their spectrum from ^^Be breakup with s, p and 
d components. 

On the other hand the ^^Be structure is the crucial input of any three- 
body model describing ^^Be as two neutrons outside an inert ^^Be core 
[101 mi US]- One then needs to know the n-^^Be interaction. Assuming 
in ^^Be a normal order of shells and an unbound s-ground state of ^^Be, the 
d5/2 resonance, experimentally at 2 MeV, has to be lowered in order to get 
the experimental two-neutron separation energy in ^^Be [HI B2] . An inver- 
sion of 2s-lpi/2 shells, similar to the inversion in ^^Be and ^°Li, was shown to 
solve this discrepancy [51 US] . This inversion is due in part to the coupling of 
the neutron with core collective excited states, inducing in the wave function 
a component with one neutron coupled to the 2"*" state of the core. Descouve- 
mont, made a GCM calculation [^1^ . He expands the ^^Be and ^^Be wave 
functions as a superposition of one and two neutrons plus a core of ^^Be, thus 
having contributions of components on excited states of the core. He gets 
an agreement with experimental values of the binding energy of ^"^Be and 
of the d5/2 resonance energy. His model gives a s-ground state in ^^Be very 
weakly bound (100 keV). However due to the uncertainties inherent to any 
model, this result may not be inconsistent with the experimental evidence 
for a weakly unbound ^'^Be. These results were confirmed later on by Baye 
and collaborators in a Lagrange-mesh calculation [3S1I1Z]- Recently Tarutina 
et al. [48j studied ^^Be and ^^Be as one and two neutrons outside a deformed 
core using a hyperspherical harmonics expansion method. They found that 
both ^"^Be and ^^Be (with an unbound s-ground state) are well described 
with a large and positive quadrupole deformation. Therefore to disentangle 
the various theoretical descriptions of ^^Be and ^'^Be one needs more precise 
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experimental information on the structure of both and in particular on their 
spectroscopic factors. 

4.4 Structure of ^^Be and ^^B 

These uncertainties in the interpretation of experimental results as compared 
to structure calculations were at the origin of our motivations to try to un- 
derstand whether the neutron-^^Be relative energy spectra obtained from 
fragmentation of ^^Be or ^^B would show differences predictable in a theo- 
retical model. It is likely that if differences will be found in the experimental 
results with ^^B and ^"^Be beams they could be due to an interplay between 
structure and reaction effects. 

The ground state of ^^Be has spin = 0^. In a simple model assuming 
two neutrons added to a ^^Be core in its ground state the wave function is: 

I'^Be >= [6i(2si/2)' + 62(1^1/2)' + &3(14/2)'] ® r'5e, 0+ > (35) 

Then the bound neutron can be in a 2s, lpi/2 or ld5/2 state. However, as 
it has been discussed in the previous section, the situation is much more 
complicated [10]- [SI and in particular the calculations of Ref. [IB] show 
that there is a large component (2si/2, 1(^5/2) ® P-Be, 2+ > with the core in 
its low energy 2"*" state which can modify the neutron distribution. 

The ground state of ^^B has spin = 2^. In a model where it is described 
as a neutron-proton pair added to a ^^Be core in its O"*" ground state with the 
proton in the IP3/2 shell, its wave function may be written as: 

\^^B >= [ai(p3/2, 2^1/2) + 02(^3/2, 4/2)] ® r'5e, 0+ > (36) 

The present experimental information [49j on ^^B is that the neutron is 
in a state combination of s and d-components with weights 66% and 30% 
respectively, while shell model calculations show a similar mixture and no 
component with an excited state of the core. There are two possibilities for 
the reaction mechanism. One is that a proton is knocked out in the reaction 
with the target. The remaining ^'^Be would be left in an unbound s-state with 
probability |aip, in a d5/2-state with probability |a2p. These unbound states 
would decay showing the s-wave threshold and d-wave resonance effects. As 
mentioned in the introduction, the second possibility is that the neutron is 
knocked out first due to its small separation energy and that the proton is 
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stripped from the remaining ^^B. We show in Sec. 5 that this can also lead 
to resonance-like effects in the cross section. 

4.5 One neutron average potential 

The link between reaction theory and structure model is made by the neutron- 
core potential determining the S-matrix in Eq. flT^ . Then if the theory fits 
the position and shape of the continuum n-nucleus energy distribution, ob- 
tained for example by a coincidence measurement between the neutron and 
the core, the parameters of a model potential can be deduced. Our initial 
bound states are obtained in a simple Woods-Saxon with R= roA^/^ 

WO = rT^(^ - 1^;;^ j -^^l + eir-R)/ay^-^ (37) 

The depth is adjusted to fit the binding energies given in Table |2] and FigJTOl 
Other parameters are also given in Table El 

Table 2: Asymptotic normalization constants Ci(fm~"'^/^) for the initial state wave func- 
tions of the bound neutron. Potential parameters are: Vq fitted to give the two ener- 
gies -0.97 MeV and -1.85 MeV, which are the known neutron binding energies in ^''B 
and in ^'^Be respectively [SD]. The other potential parameters are ro=1.27fm, a=0.75fm, 
Vso=5.25MeV. 



£.(MeV) 


-0.97 


-1.85 


li^ ji 


C.(fm- 


-1/2) 


1/2 


1.31 


1.99 


1 1/2 


0.55 


0.88 


2 5/2 


0.17 


0.34 



To describe the valence neutron in the ^^Be continuum we assume that 
the single neutron hamiltonian with respect to ^^Be has the form 

h = t + U + iW (38) 

where t is the kinetic energy and 

U{r) = Vws + SV (39) 
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is the real part of the neutron-core interaction. For the time being the imagi- 
nary part is taken equal to zero. Yws is again a Woods-Saxon potential plus 
spin-orbit whose parameters are given in Table [3l and 5V is a correction |32] : 

6V{r) = 16ae2("-^)/7(l + e^"^^^)/")^ (40) 

which originates from particle-vibration couplings. They are important for 
low energy states but can be neglected at higher energies. The above form is 
suggested by a calculation of such couplings using Bohr and Mottelson collec- 
tive model of the transition amplitudes between zero and one phonon states. 
Therefore our structure model is not a simple single-particle in a potential 
model but contains in it the full complexity of single-particle vs. collective 
couplings. However the fact that such a complexity can be put in a form like 
Eq. (l40!) is an added value to our approach. If simple fittings of experimen- 
tal data will be obtained, then the parameters of a semi-phenomenological 
potential can be deduced and linked to a more microscopic model. A more 
realistic treatment would require the description of both bound and unbound 
states in a three-body model such as in Refs. [51] and ^2] . 

Table 3: Woods-Saxon and spin-orbit potential parameters for the continuum final states. 



Vo 


ro 


a 




^SO 


(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


-39.8 


1.27 


0.75 


6.9 


0.75 



TableHlgives the energies and widths of the lpi/2 and ld5/2 states, chosen 
according to Ref.|l2] with different values of the strength a. The widths are 
obtained from the phase shift variation near resonance energy, according to 
d5j/def\e^^^ = 2/Tj, once that the resonance energy is fixed [53] • We stress 
here that while the position of our 6.5/2 resonance agrees with the experimen- 
tal evidences discussed in Sec. 4.3, the position of our lpi/2 resonance is only 
an hypothesis [121 [5i] . 

5 Results 

Results obtained with the model outlined in Sees. 2 and 4 will now be 
discussed. We describe the reaction corresponding to a neutron initially 
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Table 4: Energies and widths of unbound p- and d-states in ^'^Be and corresponding 
strength parameters for the SY potential. 





^res 

(MeV) 


(MeV) 


a 

(MeV) 


lPl/2 


0.67 


0.28 


8.34 


ld5/2 


2.0 


0.40 


-2.36 



bound in Be or B which is then excited into an unbound state of Be, 
assuming that another nucleon has been emitted and stripped by the target, 
thus not detected in coincidence with the core. The sudden approximation 
studied in Ref.[3], similar to our q=0 case, was found to be excellent for 
energy distributions like those discussed in this work. 

One of the goals of the present calculations, as far as the reaction model 
is concerned, is to understand the incident energy dependence of the breakup 
cross-section and the dependence on the neutron initial binding energy. Re- 
lated to this is the investigation of the validity of the sudden approximation 
and the accuracy necessary in calculating the phase shifts. Finally we wish to 
understand how the presence of p- and d-wave resonances, besides a threshold 
s-state in the final state, can affect the results. 

As a preliminary to a future, more accurate, study of the breakup of ^^Be 
and ^^B in a fully time dependent method, we consider the knockout of a 
single neutron from a bound state in a potential, similarly to the previous 
calculation for ^^Be. The calculations in the present section are made with 
different potentials for the initial and final state. The initial state is bound 
with a separation energy in the range 1-2 MeV and the final state is unbound. 
The continuum energies can be adjusted by varying the parameter a in the 
potential. By changing the strength a of the 6V potential in Eq. lHOj) we 
will make also the 2s-state just bound near threshold and see what would 
be the effect on the continuum spectrum. The delta-interaction strength 
used in Eq.([H]), is f 2 =-8625 MeV fm^. It has been obtained by imposing 
that this interaction gives the same volume integral as a n-^^C Woods-Saxon 
potential of strength -50.5MeV, radius 2.9 fm and diffuseness 0.75 fm. With 
a diffuseness of 0.5 fm one would obtain f2=-6717 MeV fm'^. The value 
f2=-7481 MeV fm^ would be obtained from a Woods-Saxon with the same 
geometry and a depth fitted to give the experimental neutron separation 
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sl/2tosl/2 time dependent 

s 1/2 to s 1/2 sudden 




0.5 1 1.5 2 2.5 3 

Ej (MeV) 



Figure 4: Population of resonances in the n-^^Be relative energy spectrum. 
Comparison of sudden (dashed line) and non-sudden (solid line) results for 
an s — > s transition with peak at 0.25 MeV, p^ p with peak at 0.5 MeV and 
d— > d transition with peak around 2 MeV. See text for details. 

energy in ^'^C. 

First we study the differences between results from Eqs. (!25l) and fl3T]) . the 
sudden approximation limit q=0 of those equations and Eq. (!30|) . In Figj4]we 
show the calculation of the differential probability for the transition from a 
bound s-state to an unbound s-state, a bound p-state to an unbound p-state 
and a bound d-state to an unbound d-state. The initial state parameters 
are given in Table [2j The full lines correspond to the case ej=-1.847 MeV, 
Ej„c=70 A. MeV corresponding to v=11.35 (fmxlO^^ sec~^) and q in the 
range (0.025 0.065) fm-\ The dashed lines give the q=0 calculations of 
Eqs.f l25p and ( 13T|) . There is a small difference in the absolute value of the 
probability and the sudden calculation results have slightly different widths 
in the s to s case. In the other two cases the differences are noticeable. 
We have considered only the three transitions keeping ji = jf as it would 
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a =-26.1 fm 

a =-6.6 fm 

— a=-3.5fm 
a^=-0.8fm 

- - a =4.5 fm 
-.- a=7.1fm 
a=22.4fm 
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Figure 5: Comparison of results obtained considering a final s-state for the n-^^Be relative 
energy spectrum with positive and negative scattering lengths. Scattering lengths are given 
on the figure and their corresponding SY potential strengths in Table [5] 

happen in an extreme sudden transition. Then we studied the effect of the 
extra phase in Eqs. (125|) and (!30l) on the position of the resonance peaks, and 
we found that the shift is neghgible and would not be noticeable for the range 
of neutron-core energies discussed here. Also we have calculated for several 
velocities ranging from 10 to 23 (fmxlO^^ sec~^) and found no noticeable 
differences. On the other hand changing the initial binding energy from - 
0.97 MeV to -1.85 MeV gives a widening of the distribution and a slight 
shift of the peak value, as shown in Fig. [TOl These two energies are the 
known neutron binding energies in ^^B and in ^^Be respectively [SO]- In this 
work we have kept the initial separation energy fixed at the value 1.85 MeV, 
unless otherwise stated. Our conclusion is that for fragmentation reactions 
such as those discussed here, the sudden approximation q=0 in Eqs. (125!) and 
(pTj) . is well justified from the point of view of the independence from the 
beam velocity. On the other hand the value of the Sj — Si difference has an 
important effect on the results when the final energy increases and for states 
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Figure 6: Comparison of time dependent calculation for an s to s transition with results 
of Eq.(28) using the same optical model phase shifts corresponding to scattering lengths 
as indicated. For each case we give in the legenda the normalization factor between the 
two calculations. 

with If > 0. 

The first peak shown by the experimental spectra of ^^Be needs to be 
studied in great detail as it would help determining the ground state angular 
momentum and parity. In particular, if it is due to an s-state its characteris- 
tics will depend on its closeness to threshold. Therefore we study now such 
a point. 

5.1 Low energy s-states 

Using the effective range formula [53] 

k cot ^0 = — - + -rek'^, (41) 
as 2 

in the case of a bound state of small binding energy 7 — > one has 
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Ej (MeV) 



— Eq.(28) with optical model phase shift a =-3.5 fm 

— Eq.(28) with 5=Arctan(-a^ k) 

— - Eq.(28) with phase shift from Eq.(39) 




Ej (MeV) 



Figure 7: Results of Eq.(28) using three different prescriptions for the phase shifts 
corresponding to scattering lengths as indicated. 



~- = -\l\ + ^r,^\ (42) 

Qjg / 

Equation (H2l) has been used [55] also to define an energy of unbound s-states 
of near zero energy. Such a procedure is rehable when 7/? is smaller than 
about 0.5, where R is the radius of the potential. Therefore the effective 
range formula is accurate only for very low energies. Thus we have fitted 
the behavior of our s-state phase shifts from the optical model calculation, 
to Eqs. (HT]) and ( H2i) and the corresponding parameters are given in Table [5l 
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In the case of a bound state, the effective range values can also consistently 
be obtained from Eq. 

Table 5: Strengths of the s-state i5V potential in Eq. ipD]) and corresponding scattering 
lengths, effective range parameter and energy parameter e. The strength of the central 
Woods-Saxon part is Vo=-39.8 MeV in all cases (cf. Table [3]). 



a 




re 




(MeV) 


(fm) 


(fm) 


(MeV) 


8.0 


-0.8 


117.0 




4.0 


-3.5 


17.9 




2.0 


-6.6 


11.8 




-1.0 


-26.1 


7.58 




-5.0 


22.4 


5.9 


0.06 


-15.0 


7.1 


3.8 


1.34 


-35.0 


4.5 


2.7 


6.49 



Figure [5] shows the influence of the phases 5 and v (cf. Eq.( IT7|) ) on the 
breakup cross sections. The results shown correspond to final s-states with 
positive and negative scattering lengths. Several cases are considered and the 
corresponding potentials, scattering lengths, and effective ranges are given in 
Table [51 The scattering length values were obtained from the phase shifts as 
as = — lim , and also cross-checked by solving the Schrodinger equation 

at zero energy. The bound state energies in the last column were obtained 
from the solution of the Schrodinger equation in each given potential. Notice 
that the breakup cross sections for potentials with negative scattering lengths 
are much larger than those with positive scattering lengths. This effect is 
mainly due to the influence of the phase u. 

The effective phase in Eq.( ll7p is 5 = 6 + u and there are interference 
effects which are especially important for an s-state final state. When k is 
small S ~ —ktts and it depends on the sign of the scattering length, while 
z/ ^ kipc -|- 1/7) and it is always positive. The part of the cross section 
Eq. (l3Tl) with the probability Eq. (l25|l for = and small /c, depends on the 
relative sign of as and v as 

^ (x\5 + ^\-kas + v\^. (43) 

When a^ < the cross section is increased relative to the value at = 0, 
while for a^ > it is reduced. This effect is seen clearly in FigO These 
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interference effects can also explain why the cross section for as=4.8 fm is 
larger than for as=7.1 fm. Also, the decrease in the positive values of a^ 
shown in Table [5] corresponds to an increase in the depth of the potential 
6V. Fig. [5] shows that the cross section increases for the smallest values of 
a^. As the attraction becomes ever stronger the scattering length changes 
sign and the cross section becomes larger. This effect is due to the influence 
of a continuum s-state coming close to threshold. When the final potential 
has a very weakly bound 2s-state with as=22.4 fm one sees a very narrow 
peak close to threshold (thick solid line) while for as= 4.5 fm, corresponding 
to a more strongly bound state, no peak at all, rather a kind of flat bump 
( double- dot dashed line ) . 
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Figure 8: Results obtained including the s, p and d states. Each curve corresponds to 
just one transition as indicated. The solid curve is the sum of all transitions from the 
s-bound state. To make them visible some curves have been multiplied by a factor of five 
as indicated in the legenda. 

The calculations presented in Fig. 5 of Ref . [7] , did not include the extra 
phase u because the final state interaction with the core of origin was ne- 
glected while the final state interaction with respect to the other nucleus was 
included. In that case the neutron behaved as a free particle in the scattering 
on the other nucleus. Breakup cross sections depended quite strongly on the 
magnitude of the scattering length but had a weak dependence on its sign. 
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The reaction mechanism discussed in this paper is instead an inelastic-type 
of excitation in which the final state interaction is with the core of the projec- 
tile and therefore the present formalism shows that the S-matrix in Eg. (1191) 
as well as in Eq. fl30|) is effectively off-the-energy-shell. In the s ^ s case we 
show also in Figj6l a comparison between the results just discussed and those 
obtained using the optical model phase shifts in Eq. fl5(I]) and whose absolute 
values have been normalized to our peaks. As anticipated in Sec. 2.3 and 2.4, 
the curves from the time dependent model show a faster decrease towards 
high energies. This is because the form factors in Eq. (l25|) decrease rapidly 
and they have large values only for impact parameters close to the strong 
absorption radius. 




Figure 9: Check of the dependence from the initial state angular momentum. Full curve: 
sum of transitions from s-initial state. Dashed and dotdashed lines: sum of transitions 
from p and d-initial states respectively. To make them visible these curves have been 
multiplied by a factor of four and eight respectively. 



Since Eq.f l30p is quite simple to implement in an fit to experimental spec- 
tra, we have also studied its sensitivity to the choice of the phase shift. FigH 
shows, for two scattering lengths, results obtained using optical model phase 
shifts, the second order effective range approximation Eq. fHTj) with values 
from Table [5l and the first order phase shift 6 ~ indicated. As 

expected, the latter approximation is reliable only for extremely small values 
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of the final energy. The second order, effective range parametrization works 
much better, in particular as the scattering length increases. 

5.2 p and d- resonances 

Fig. [S] shows results obtained considering three different possibilities for the 
initial state: the s, p and d orbitals. For each initial state a unit spectroscopic 
factor is assumed. 




Figure 10: Check of the dependence from the initial binding energy of the sum of 
transitions from s-initial state. FuU curve: Ei— -0.97 MeV as in ^''B; dashed curve: ei=- 
1.85 MeV as in ^''Bc [50 . Dotdashed hue: sum of transitions from s, p and d- initial states 
including spectroscopic factors of 0.66, 0.04. 0.30 respectively as in ^^B [49] with ei= -0.97 
MeV. 



We show the results of the transition bound to unbound from each initial 
state to each possible unbound state as indicated. Available experimental 
spectra [23l [26] show that the next group of resonances is located around 4-5 
MeV. For this reason higher partial waves have not been included. We have 
checked that if the d3/2 and P3/2 states are calculated in the same potentials 
as those of the d5/2 and pi/2, they would give a noticeable, non resonant, 
contribution only for a transition from the s initial state. This contribution, 
not shown in the figures, would constitute a smooth background. The overall 
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effect would be an increase the cross section value of about 10% around the 
0.5 MeV peak and of about 25% around 2 MeV in FigJHl The spectrum for 
Coulomb breakup has also been calculated and found negligible, compared to 
the other transitions, for the separation energies in ^^Be and ^^B. Therefore 
it is not included in the figures. 

The p and the d-resonance peaks are clearly seen because of the effect of 
the angular momentum enhancement factor in Eq. flS^ . As indicated some 
transition strengths have been multiplied by a factor five to make them visi- 
ble. It is clear from this figure that the dominant components in the neutron 
knockout spectrum from ^^Be and ^^B come from the s-wave component in 
the ground states of those nuclei. Therefore the full curve is the sum of the 
contributions from the initial s-state alone with unit spectroscopic factor. 
There can be large peaks due to transitions to the pi/2 and 6.5/2 final state, 
provided they are centered around the energies we have choosen. The results 
of Fig. [9] are shown to check the dependence of the transition probability on 
the initial state angular momentum. The full curve is the sum of transitions 
from s-initial state. Dashed and dotdashed lines are the sum of transitions 
from p and d-initial states respectively. Since the transitions from p and d 
orbitals are negligible, then components in the wave functions of ^^Be or ^^B 
with a neutron in such angular momentum states will not play much role in 
the reaction. Thus it is unlikely that any difference in the neutron breakup 
spectrum due to different mixtures of these configurations in the two parent 
nuclei will be seen. 

To clarify further the latter point, the sum of all transitions from the s- 
bound state, is shown again in Fig. [10] for two initial binding energies (solid 
and dashed lines as indicated). Very small changes in the initial binding 
energy do not give appreciable differences in the final continuum spectra, in 
particular in the positions of the peaks. They however give differences in the 
absolute cross section value and in the relative height of the s- and d-state 
peak. The dotdashed line is the result obtained using the neutron initial 
binding energy in ^^B, which is -0.969 MeV, and summing transitions from 
s, p and d-initial states including experimental spectroscopic factors [IH] of 
0.66, 0.04. 0.30 respectively. According to the simple model for two nucleon 
breakup presented in Sec. 3, the presence of a lp3/2 proton in the initial state 
would give the same contribution to the breakup of each neutron component 
and therefore the shape of the spectra should not be modified. On the other 
hand in the case of ^"^Be the two neutrons are in the same state for each 
component of the initial wave function. Therefore since the p and d wave 



31 



functions have less pronounced tails, the last integral of Eg. (1551) . which gives 
the stripping probability of one of the two neutrons, will naturally diminish 
the absolute value of the p and d resonances peaks with respect to the final 
s-state peak. 

We notice also that there is a shift between the peaks of the energy spectra 
and the resonance energies obtained from the phase shifts and given in Table 
4. The shift is due to a combined effect of the 1/k factor in Eq. fE^ . of the 
matching between initial separation energy, peak energy of the resonance 
and relative beam energy per nucleon and only in a small measure to the 
presence of the extra phase v. The matching effect is manly contained in 
the slope of the form factor, which to first order is equal to the decay length 
of the initial state wave function (cf. Eq. (l2^ ). In the ^^Be case it is less 
noticeable because the initial separation energy is very small. The resonance 
energy was 1.27 MeV while the peak is at 1.25 MeV. For ^^Be the resonance 
energies are 0.67 MeV (p-state) and 2 MeV (d-state) while the peaks are at 
0.5 MeV and 1.8 MeV respectively. One notices also that the shift increases 
increasing the resonance energy. 

Therefore two differences might be expected in the experimental energy 
spectra of ^^Be when produced by ^^Be rather than ^^B projectile fragmen- 
tation: a first peak at energy smaller than 0.5 MeV due to the s-continuum 
state. The s-continuum peak below 0.5MeV, as given by our calculations, can 
be seen better in Fig. 8 and it is due to the s to s transition. If the s initial 
component would contribute in reality more than the p and d components 
(we have taken them with equal weights) such a peak could be more evident 
in the data. Then two more diffuse bumps if the projectile is ^^Be. Two well 
definite peaks of not too different strength, one centered at 0.5 MeV due to 
the pi/2 resonance and another around 2 MeV due to the resonance if 
the projectile is ^"^B. Such an hypothesis agrees also with the conclusion of 
Ref. [14j of a s-state very close to threshold. The three different experiments 
[5l [HI [15] would therefore be complementary and allow to determine the 
characteristics of ^^Be low energy continuum. 

To give another example of a possible comparison with available data, we 
show in Fig. [TT] the experimental points from H. Simon et al. [SS] for the 
reaction ^^Be+^^C n+^^Be+X at 250 A. MeV. The normalization factor of 
the data to mb/MeV is 0.843. The solid line gives the sum of all transitions 
from the s initial state with £:/=-1.85 MeV (solid line), as in previous figures, 
renormalized with a factor 2.4. The dashed line is the folding of the calculated 
spectrum with the experimental resolution curve. Therefore the calculation 
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Figure 11: Sum of all transitions from the s initial state with £^=-1.85 MeV (solid line). 
Experimental points from H. Simon et al. for the reaction ^"^Be+^^C n+^^Be+X at 
250 A. MeV. Dashed line is the folding of the calculated spectrum with the experimental 
resolution curve. 

underestimate the absolute experimental cross section by a factor of 2. In 
view of the incertitude in the strength of our n-target (5-potential and on 
the initial state spectroscopic factor which has been taken as unit, we can 
consider our absolute cross sections quite reasonable. 

Finally we wanted to address the issue of possible core excitation effects 
which in Ref. [58] have been shown to be of fundamental importance for repro- 
ducing simultaneously the ^^Be and ^^Be characteristics. Those effects can 
be modeled in the present approach by considering a small imaginary part 
in the neutron-core final optical potential (cf. Eq. fISS]) ). This is a standard 
procedure for continuum states where most often the potential is also energy 
dependent |57]. A surface potential of Woods-Saxon derivative form has been 
taken with strengths Wq equal to -0.5, -1.0 and -1.5 MeV for the d-state only. 
The results are shown in Fig. [121 The effect of the imaginary potential is 
to wash out the d-resonance peak. Several structure models predict indeed 
the d5/2 resonance coupled to an excited ^^Be core. We have found the same 
smoothing off effect if the s-state is calculated including an imaginary poten- 
tial. It seems therefore that the spectrum of unbound nuclei could partially 
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reflect the structure of tlie bound parent nucleus and that reaction mecha- 
nism models used to extract structure information should carefully include 
the effects discussed above. The model presented here seems to be quite 
promising in this respect. 




3 



(MeV) 

Figure 12: Sum of all transitions from the s initial state with 85 MeV including 

core excitation via an imaginary part of the optical potential for the d-resonance only. 



6 Conclusions and outlook 

In this paper we have presented a model to study one neutron excitations 
from a bound initial state to an unbound resonant state in the neutron-core 
low energy continuum. This is the process by which unbound nuclei are 
created and studied via projectile fragmentation experiments [5]- pn [56] . 

The model is based on a time dependent perturbation theory amplitude 
and the final state is described by an optical model S-matrix. It can be 
considered an evolution with respects to sudden and/or R- matrix theory 
models. The advantages are that the model can be applied to fragmentation 
from deeply bound states and to resonant and non resonant, large energy, 
continuum final states. Also core excitation effects can be modeled by an 
imaginary part of the neutron-core optical potential. 
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A related approach has been developed some time ago and applied to the 
treatment of transfer to the continuum in which, following the interaction 
between two passing-by nuclei, a neutron from one of them goes from a 
bound to an unbound state with final state interaction with the other nucleus. 
Comparison of the present formalism to the transfer to the continuum model 
shows that in principle projectile fragmentation does not reflect directly the 
properties of the neutron-core resonances because the reaction mechanism 
induces an extra phase with respect to the free particle neutron-core phase 
shift. It means that the measurements would probe an off-the-energy-shell 
S-Matrix. The distortion effects seem however small and negligible for the 
cases discussed in this work. 

One neutron breakup can be studied in this way but also one step of two 
neutron breakup of a borromean nucleus. In this paper we have presented 
some applications to both cases to study the properties of ^^Be continuum 
and of ^^Be. Our results are in agreement with the conclusions of Ref . [U [39] 
for ^^Be. Due to the structure inputs we use, in particular the position of the 
Pi/2 resonance, the ^^Be continuum spectrum obtained from fragmentation of 
^■^B or ^^Be shows essentially the effect of the continuum p and d-resonances. 
The s-state although present in the calculations almost disappears inside the 
tail of the p-state but it would still determine the ground state spin and 
parity of ^^Be. Obviously we cannot be conclusive on the structure of ^^Be 
because at the moment we have not attempted to fit experimental data but 
simply to develop a good reaction model. Furthermore our structure inputs, 
although reasonable, are extremely simple compared to the complexity of the 
nucleus under study. However preliminary comparisons seem to indicate the 
reliability of our model. 

We have also shown that the excitation energy spectra of an unbound 
nucleus might reflect the structure of the parent nucleus from whose frag- 
mentation they are obtained. In particular, in the case of ^"^Be fragmentation, 
the initial state spectroscopic factors are not known experimentally, and the 
information from structure calculations indicate an important configuration 
mixing with components coupled to an excited ^^Be core. Thus the analysis 
of such spectra is expected to be even more complicated. 

However from the point of view of reaction theory, the results obtained 
here seem promising and we hope to use such a procedure to implement a 
fit of experimental data on unbound nuclei. At the same time we intend 
to construct an accurate, second order, fully time dependent theory of two 
neutron breakup, incorporating properly the time ordering between the two 
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neutrons. In this way we hope to clarify the question of sequential versus 
simultaneous mechanism implicit in the formation of neutron-core resonance 
states in reactions like ^^Li + X ->^°Li ->^Li+2n [5j or ^^Be + X ^^^Be*+n 
+ X^^^Be+2n + X, or ^^B + X -^^^ Be*+p + X^^^Be+n +p + X. 
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A Modifications to the (5- interact ion 



Figure 13: Graph of variables used in the calculation of appendix A. 




The purpose of this appendix is to justify the use of a (5-interaction as an 
approximation for the finite range n-target interaction and to derive Eq. (|TOl) . 
We then calculate 

J= dxdydzdz' — e-"'^'-'^V2(x-bc,y,z'). (44) 

J-oo 

If 7 is large the integral is concentrated near the surface of V2(r). To simplify 
the discussion put q = 0. Also, as in Sec. 2.3 put 

g-(7-jfc)r ^ ^-ax ^-a{y'^+z'^)/2x 



^-ax^-azy2Ro^ (45) 
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where a = 'j — ik and we used the fact that the gaussian term e""*^^ """^ -'''^^ 
gives the largest contribution at Rq = be — Rt which is the position of the 
surface of V2(r). To simphfy further the calculation we have neglected the 
y-dependence in Eq. fHSj) but kept the z-dependence so that the integral J 
will converge. Indicate x = x — Rq = x — (be — Rt)- Take V2(r) to be a 
square well potential of depth Vq and radius Rt- Then 

/oo 
1^2(3; - be, y, z') dydz' = nVof3^ = 2WqRtx (46) 
-00 

where the upper limit of the two dimensional integral is given by + (x — 
Rt? = Rt and /J^ ^ 2xRt- 
So that 

J = / cixe-""x / dze~^ 



hi 



bl 



'27r(6e 


-Rt) 


7 - 


- ik 


'27r(6e 





27r:^e-"(''^-^^) (47) 



V '~) — ik 

The ratio of the integral J to the integral I of Eq. fl27|) is 



1 



/ \ b, J 2R^{-f-ikf' 

Where we imposed that the strength of the 5-interaction be equal to the 
volume integral of the square well potential V2 = ^tiVoRt- Thus to represent 
the effect of a finite range potential by a 5- interaction when 7RT >> 1, 
replace 

• (1) be ^ be — Rt i-e. the interaction is at the surface of the target. 

• (2) Multiply the strength of the interaction by | (y-ik)^ ■ 

This factor is less than one. The change (1) increases the breakup integral, 
the factor (2) decreases it. 
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B Including spin 

To include spin variables in the initial and final states is mainly an angular 
momentum coupling problem. The angle-spin wave function of the initial 
and final states are 

"^iUi, li, 0,(f)) = {jini\hmi^ai)Yi.^X^, (fyXaXp) (49) 

'^fijf,nf,lf,e,(f)) = (3fnf\lfmf^af)Yi^.rnf{0,(l))Xaf{p)- (50) 

We choose the quantization axis along the y-direction, such that = 0. Then 
after integration over p, the angle spin part of the overlap Eq.(l9]) is: 

1 1 

V(jfnf,j,ni) = Y (i/'^/l^/"^/^^)(ii^il^i"^i^^)^i*rn^(^'0)^'.m.(^'0) 

= E (-l)™^(j/ri/|//-m^^a)0-n,|/,m,^a)y,,„,(e,0)y^,„^(e,0), (51) 

mfmiu 

where we have put = {—l)"^Yi_m- Next we use the relation for coupling 
two spherical harmonics of the same argument and introduce the notation 
I = V2z + 1 

Yi^^XO,0)Yi..^je,0) = Y{LM\kmilfmf){LO\kOlfO)^^YLM{0,0). (52) 

LM v47rL 

Substituting into the relation for V{jfnf,jini) there is a sum of products of 
three Clebsch-Gordan coefficients which reduces to a product of a Clebsch- 
Gordan coefficient and a 6-j symbol. Collecting together the terms evaluated 
above we get: 

'DUfnf,j,n,) = (-)fY{LM\j,n,jj-nf){LO\mfO)^-^^Y,M{e,0) ( ^ 

LM V^TtL [ Jf Ji 

With the phase f = nj + If — j j 

In this scheme the integral Eq. (l26|) 

rf^e^''^^S/^r(^7r)lz„™,(^,0)A;-/iS:)(A;r)l^^,^^(e,0) (53) 

~oo Z 
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is substituted by a new integral I^m defined as 

Ilm = r dze''^H'^^hl^\i^r)k'-h\-\kr)YL,Mi(^, 0). (54) 

Summing over nf and averaging over and using the orthogonality of the 
{LM\iinii fHf) coefficients while calculating |A(jjn/, as in Eq. (l22l) we 

find that Eq. (l2H]) can be replaced by: 



(55) 

where 

r(i ■ 1 ■ n (2j/ + l)(2/. + !)(% + !) f /, If lV 

A sum rule for 6 — j symbols gives 

i:C(l.M.L)-ilg^y (57) 

On the other hand, if in Eq. (!25i) and ( l26l) or ( l53l) we use the angular 
momentum coupling formula for spherical harmonics 

^ {LM\km^lfmJ)Yl^„^X0,O)Yl^„,^{e,O) = (L0|/,0//0)4^>LAf(^, 0), 

niimf V 47rL 

(58) 

then the relation between Im mt and Itm is 

Ikmafm, = Y.i-'^r'{LM\lf-mfkm,){LO\hOlfO)^^lLM (59) 
LM V47rL 

and Eq. (l25i) is replaced by 

which could also be obtained using Eq. flFTl) in fl5Sl) . 
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